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Abstract 

During B cell development, the precursor B cell receptor (pre-BCR) checkpoint is thought to increase immunoglobulin k 
light chain [IgK) locus accessibility to the V(D)J recombinase. Accordingly, pre-B cells lacking the pre-BCR signaling 
molecules Btk or Slp65 showed reduced germline transcription. To investigate whether pre-BCR signaling modulates Vk 
accessibility through enhancer-mediated IgK locus topology, we performed chromosome conformation capture and 
sequencing analyses. These revealed that already in pro-B cells the k enhancers robustly interact with the ~3.2 Mb Vk 
region and its flanking sequences. Analyses in wild-type, Btk, and Slp65 single- and double-deficient pre-B cells 
demonstrated that pre-BCR signaling reduces interactions of both enhancers with IgK locus flanking sequences and 
increases interactions of the 3'k enhancer with genes. Remarkably, pre-BCR signaling does not significantly affect 
interactions between the intronic enhancer and genes, which are already robust in pro-B cells. Both enhancers interact 
most frequently with highly used Vk genes, which are often marked by transcription factor E2a. We conclude that the k 
enhancers interact with the Vk region already in pro-B cells and that pre-BCR signaling induces accessibility through a 
functional redistribution of long-range chromatin interactions within the Vk region, whereby the two enhancers play distinct 
roles. 
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Introduction 

B lymphocyte development is characterized by stepwise 
recombination of immunoglobulin (Ig), variable (V), diversity 
(D), and joining (J) genes, whereby in pro-B cells the Ig heavy (H) 
chain locus rearranges before the IgK or IgX light (L) chain loci 
[1,2]. Productive IgH chain rearrangement is monitored by 
deposition of the IgH \x chain protein on the cell surface, together 
with the preexisting surrogate light chain (SLC) proteins X5 and 
VpreB, as the pre-B cell receptor (pre-BCR) complex [3]. Pre- 
BCR expression serves as a checkpoint that monitors for 
functional IgH chain rearrangement, triggers proliferative expan- 
sion, and induces developmental progression of large cycling into 
small resting Ig ll"*" pre-B cells in which the recombination 
machinery is reactivated for rearrangement of the IgK or IgA L 
chain loci [3,4]. 



During the V(D]J recombination process, the spatial organiza- 
tion of large antigen receptor loci is actively remodeled [5]. 
Overall locus contraction is achieved through long-range chro- 
matin interactions between proximal and distal regions within 
these loci. This process brings distal V genes in close proximity to 
(D]J regions, to which Rag (recombination activating gene) protein 
binding occurs [6] and the nearby regulatory elements that are 
required for topological organization and recombination [5,7,8]. 
The recombination-associated changes in locus topology thereby 
provide equal opportunities for individual V genes to be 
recombined to a (D)J segment. Accessibility and recombination 
of antigen receptor loci are controlled by many DNA-binding 
factors that interact with local m-regulatory elements, such as 
promoters, enhancers, or silencers [7-9]. The long-range chro- 
matin interactions involved in this process are thought to be 
crucial for the regulation of V(D)J recombination and orchestrate 
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Author Summary 

B lymphocyte development involves the generation of a 
functional antigen receptor, comprising two heavy chains 
and two light chains arranged in a characteristic "Y" shape. 
To do this, the receptor genes must first be assembled by 
ordered genomic recombination events, starting with the 
immunoglobulin heavy chain (IgH) gene segments. On 
successful rearrangement, the resulting IgH |i protein is 
presented on the cell surface as part of a preliminary 
version of the B cell receptor — the "pre-BCR." Pre-BCR 
signaling then redirects recombination activity to the 
immunoglobulin k light chain gene. The activity of two 
regulatory k enhancer elements is known to be crucial for 
opening up the gene, but it remains largely unknown how 
the hundred or so Variable (V) segments in the k locus gain 
access to the recombination system. Here, we studied a 
panel of pre-B cells from mice lacking specific signaling 
molecules, reflecting absent, partial, or complete pre-BCR 
signaling. We identify gene regulatory changes that are 
dependent on pre-BCR signaling and occur via long-range 
chromatin interactions between the k enhancers and the V 
segments. Surprisingly the light chain gene initially 
contracts, but the interactions then become more func- 
tionally redistributed when pre-BCR signaling occurs. 
Interestingly, we find that the two enhancers play distinct 
roles in the process of coordinating chromatin interactions 
towards the V segments. Our study combines chromatin 
conformation techniques with data on transcription factor 
binding to gain unique insights into the functional role of 
chromatin dynamics. 

changes in subnuclear relocation, germline transcription, histone 
acetylation and/ or methylation, DNA demethylation, and com- 
paction of antigen receptor loci [5,10]. 

The mouse IgK locus harbors 101 functional Vk genes and four 
functional elements and is spread over >3 Mb of genomic 
DNA [11]. Mechanisms regulating the site-specific DNA 
recombination reactions that create a diverse IgK repertoire are 
complex and involve local differences in the accessibility 
of the Vk and Jk genes to the recombinase proteins [12]. 
Developmental-stage-specific changes in gene accessibility are 
reflected by germline transcription, which precedes or accom- 
panies gene recombination [13]. In the IgK locus, germline 
transcription is initiated from promoters located upstream of Jk 
(referred to as k" transcripts) and from Vk promoters [14]. 
Deletion of the intronic enhancer (IEk), located between Jk and 
Ck, or the downstream 3'k enhancer (3'Ek), both containing 
binding sites for the E2a and Irf4/Irf8 transcription factors (TFs), 
diminishes IgK locus germline transcription and recombination 
[15-19]. On the other hand, the Sis (silencer in intervening 
sequence) element in the Vk^k region negatively regulates IgK 
rearrangement [20]. This Sis element was shown to target IgK 
alleles to centromeric heterochromatin and to associate with the 
Ikaros repressor protein that also colocalizes with centromeric 
heterochromatin. Sis contains a strong binding site for the zinc- 
finger transcription regulator CTCC-binding factor (Ctcf) 
[21,22]. Interestingly, deletion of the Sis element or conditional 
deletion of the Ctcf gene in the B cell lineage both resulted in 
reduced k'^ germline transcription and enhanced proximal Vk 
usage [21,23]. Very recently, a novel Ctcf binding element 
located directly upstream of the Sis region was shown to be 
essential for locus contraction and recombination to distal Vk 
genes [23]. In addition, the IgK repertoire is controlled by the 
polycomb group protein YYl [24]. 



Induction of IgK rearrangements requires the expression of the 
Ragl and Rag2 proteins, the attenuation of the cell cycle, and 
transcriptional activation of die IgK locus, all of which are thought 
to be crucially dependent on pre-BCR signaling [4,25]. At first, 
pre-BCR signals synergize with interleukin-7 receptor (IL-7R) 
signals to drive proliferative expansion of IgH large pre-B cells 
[4] . In these cells, transcription of the Rag genes is low and the 
Rag2 protein is unstable due to cell-cycle-dependent degradation 
[26]. Subsequently, signaling through the pre-BCR downstream 
adapter Slp65 (SH2-domain-containing leukocyte protein of 
65 kDa, also known as Blnk or Bash) switches cell fate from 
proliferation to differentiation [4]. Importandy, Slp65 (i) induces 
the TF Aiolos, which down-regulates Xb expression [27]; (ii) binds 
Jak3 and thereby interferes with IL-7R signaling [28]; and (iii) 
reduces inhibitory phosphorylation of Foxo TFs [29]. All these 
changes result in attenuation of the cell cycle and thus Rag protein 
stabilization. Moreover, Rag gene transcription is induced by Foxo 
proteins [30]. 

Although rearrangement and expression of the IgK locus can 
occur independently of IgH |i chain expression [31,32], several 
lines of evidence indicate that pre-BCR signaling is actively 
involved in inducing IgK and IgX locus accessibility and gene 
rearrangement. First, surface IgH |i chain expression correlates 
with germline transcription in the IgK locus [33]. Second, in the 
absence of Slp65, k" germline transcription is reduced [34]. 
Third, mice deficient for Bruton's tyrosine kinase (Btk), which is 
a pre-BCR downstream signaling molecule interacting with 
Slp65, show reduced IgX L chain germline transcription and 
reduced Igk usage [35]. Fourth, transgenic expression of the 
constitutively active E41K-Btk mutant in IgH |t chain negative 
pro-B cells induces premature rearrangement and protein 
expression of IgK L chain [34]. Based on fluorescence in situ 
hybridization (FISH) studies, it has been proposed that in pro-B 
cells distal Vk and Ck genes are separated by large distances 
and that the IgK locus specifically undergoes contraction in 
small pre-B and immature B cells actively undergoing Vk-Jk 
recombination [36]. However, it remains unknown how pre- 
BCR-induced signals affect the accessibility, contraction, and 
topology of the Vk region, or how they affect the long-range 
interactions of the K regulatory elements involved in organizing 
these events. 

In this study, we identified the effects of pre-BCR signaling on 
germline Vk transcription and on the expression of TFs 
implicated in the regulation of IgK gene rearrangement. We 
found that the decrease in pre-BCR signaling capacity in wild- 
type, Btk-deficient, Slp65-deficient, and Btk/Slp65 double- 
deficient pre-B cells was paralleled by a gradient of decreased 
expression of many TFs including Ikaros, Aiolos, Irf4, and (to a 
lesser extent) E2a, as well as by a decreased IgK locus accessibility 
for recombination. Several of these factors can mediate long- 
range chromatin interactions and are known to occupy K 
regulatory elements that regulate locus accessibility [37-40]. 
We therefore sought to analyze the effect of pre-BCR signaling 
on the higher order chromatin structure organized by these 
regulatory sequences at the IgK locus. To this end, we performed 
chromosome conformation capture and sequencing (3C-seq) 
analyses [41] on pro-B cells and pre-B cells from mice single or 
double deficient for Btk or Slp65 to evaluate the effects of this 
pre-BCR signaling gradient on IgK locus topology. These 3C-seq 
experiments demonstrated that already in pro-B cells the K 
enhancers robustiy interact with the ~3.2 Mb Vk region and its 
flanking sequences, and that pre-BCR signaling induces accessi- 
bility by a functional redistribution of enhancer-mediated 
chromatin interactions within the Vk region. 
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Figure 1 . Gene expression profiling strategy for the identification of genes regulated by Bt1</Slpe5-mediated pre-BCR signaling. (A) 

FACS sorting strategy for purification of pre-B cell fractions from the indicated mice on a VH81x transgenic Ragl^'^ background. Lymphocytes were 
gated on the basis of forward/side scatter and B220*CD19* pre-B cell fractions were sorted. Virtually all B220*CD19* cells were cytoplasmic \i heavy 
chain positive [34], but showed genotype-dependent levels of expression of the CD2 differentiation marker, in agreement with previous findings [34]. 
(B) DNA microarray analysis of total mRNA from FACS-purified B220*CD19* pre-B/pro-B cell fractions from the indicated mice. One-way ANOVA 
analysis (p = 0.01) identified 266 significantly different genes. IVleV hierarchical clustering of gene expression differences are represented in the 
heatmap. (C) Validation of the expression of TFs implicated in IgK gene rearrangement. Total mRNA isolated from FACS-sorted B220*CD19* pre-B/ 
pro-B cell fractions from the indicated mice was analyzed by quantitative RT-PCR for expression of TFs. Expression levels were normalized to those of 
Gapdh, whereby the values in WT pre-B cells were set to one. Bars represent mean values and error bars denote standard deviations for four 
independent mice per group. 
doi:10.1371/journal.pbio.1001791.g001 



Results 

Identification of Genes Regulated by Pre-BCR Signaling 

Whereas mice deficient for the pre-BCR signaling molecules Btk 
and Slp65 have a partial block at the pre-B cell stage [42,43], in 
Btk/Slp65 double-deficient mice, only very few pre-B cells show 
progression to the immature B cell stage characterized by functional 
IgL chain gene recombination [44]. To enable analysis of the effects 
of pre-BCR signaling on (i) the expression of genes involved in IgK 
gene rearrangement and on (ii) long-distance chromatin interactions 
in the IgK locus in pre-B cells in the absence of IgK gene 
recombination events, we bred Btk and Slp65 single- and double- 
deficient mice on the Ragl ' background. In these mice, 
progression of B cell progenitors to the pre-B cell stage was 
conferred by the transgenic, functionally rearranged VH81x IgH |i 
chain, which ensures pre-BCR expression and cellular proliferation. 



The absence of functional Ragl protein precludes IgL chain gene 
rearrangement and cells are completely arrested at the small pre-B 
cell stage (Figure lA). 

We performed genome-wide expression profiling of FACS- 
purified B220"^CD 1 9"" pre-B ceU fractions from wild-type (WT), Btk, 
and Slp65 single- and double-deficient VH81x transgenic Ragl~'~ 
mice (Figure lA). In these experiments non-VH81x transgenic 
Ragl~'~ pro-B cells served as controls. One-way ANOVA analysis 
using MeV software (/)<0.0I) [45] revealed that 266 genes were 
diflferentially expressed between the five groups of pro-B/pre-B cells 
(Figxire IB). When compared with WT VHBlx transgenic Ragl~'~ 
pre-B cells, 174 genes were up-regulated, whereby the average 
values of the fold increase were ~1.70, ~3.28, ~3.36, and ~3.47 
for Btk~'~, Slp65~'~,Btk~'~Slp65~'~ VH81x transgenic 
pre-B cells and non-VH81x transgenic Ragl~'~ pro-B cells, 
respectively (see Table SI). A similar gradient of gene expression 
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changes was appairent from the average values of the fold change for 
the 192 significantly down-regulated genes, which were ~1.65, 
~2.29, ~3.79, and ~4.15 in the four groups of pre-B/pro-B cells, 
respectively (see Table S2). In a hierarchical clustering analysis of 
the five groups of B cell precursors, the expression profiles oiBtk ' 
Slp65~'~ VH81x transgenic Ragl''' pre-B cells and non-VH81x 
transgenic Ragl ' pro-B cells were very similar (Figure IB). This 
implies that expression of the 266 genes is not substantially 
influenced by pre-BCR-mediated proliferation, which is still 
induced in pre-B cells lacking both Btk and Slp65 [44,46] but not 
in Ragl ' pro-B cells. Consistent with these findings, gene distance 
matrix analysis revealed a clear gene expression gradient among the 
five groups of pre-B/pro-B cells, in which BA '~Slp65 ' pre-B 
and Ragl~'~ pro-B cells again showed highly comparably 
expression signatures (Figure SI). 

In agreement with previous findings [,34,43,46], pre-BCR 
signaling-defective pre-B cells manifested increased expression of 
Dntt, encoding terminal deoxynucleotidyl transferase and the SLC 
components Vpre (Vprebl) and A,5 {Iglll), as well as decreased 
expression of the cell surface markers Cd2, Cd22, Cd25(IL-2R), 
and MHC class II (Table 1). Btk and Slp65 single-deficient and 
particularly double-deficient pre-B cells failed to up-regulate 
various genes known to be involved in IgL chain recombination, 
such as Ikzfi (Aiolos), Ikzfl (Ikaros), Iif4, Spib, Pou2J2 (Oct2), 
polymerase-n [47], as well as Hivepl encoding the Mbp-1 protein, 
which has been shown to bind to the K enhancers [48]. In 
addition, pre-BCR signaling influenced the expression levels of 
many other DNA-binding or modifying factors that were not 
previously associated with IgL chain recombination, including 
Lmo4, ^p710, Aridla/3a/3h, the lysine-specific demethylases Aofl 
and Phf2, Prdm2 (a H3K9 methyltransferase), the sikl gene 
encoding a histon deacetylase (HDAC) kinase, Hdac5, Hdac8, 
and the DNA repair protein gene Revl (Table S2). We did not find 
significant dilferences in the expression of several other TFs 
implicated in Ig gene recombination — for example, Obf 1 / Oca-B, 
Pax5, E2a, and IrfB (Table 1). In addition, in signaling-deficient 
pre-B cells, we found reduced transcription of genes encoding 
several signafing molecules (e.g., Rasgrpl, Rapgefll, Ralgps2, Blk, 
Traf5, Hck, Nfkhia (iKBa), Syk, Csk), cell surface markers (Cd38, 
Cd72, Cd74, Cd55, and Notch2), or genes regulating cell sur\-i\^al 
(Bnif and Bcl2ll encoding BcIxl) (Table S2). Interestingly, we 
observed concomitant up-regulation of signaling molecules that 
are also associated with the T cell receptor [Lat, Z'^pJO, and Prhq 
(PKCe); Table SI). 

Next, we used quantitative RT-PCR to confirm the observed 
differential expression of several TFs. Expression levels of these 
genes were indeed significantly reduced in a pre-BCR signaling- 
dependent manner, especially for Aiolos, Ikaros, and Irf4, with 
residual expression levels in Btk~'~ Slp65~^~ VHBlx transgenic 
Ragl~'~ pre-B cells that were —1%, —20%, and —9% of those 
observed in VVT VH81x Ragl~'~^ mice, respectively (Figure IC). 
In addition, we found moderate effects on Obfl (Oca-B) and E2a 
with residual expression levels of ~28% and ~44%, respectively. 
In chromatin immunoprecipitation (ChIP) assays, we observed in 
pre-B cells substantial binding of E2a protein to the intronic and 3' 
K enhancer regions and to the three Vk regions analyzed. Under 
conditions of reduced pre-BCR signaling activity, E2a binding to 
the enhancers was essentially maintained (3'Ek) or reduced (IEk), 
but E2a binding to the Vk regions was lost (Table S3). Consistent 
with the significant reduction of Ikaros expression in 5'^ 55 ' pre- 
B ceUs, Ikaros binding to both k enhancers and Vk regions was 
undetectable in these cells (Table S3). 

Taken together, from these findings we conclude that the 
five groups of pro-B/pre-B cells, representing a gradient of 



progressively diminished pre-BCR signaling, show in paraUel a 
gradient of diminished modulation of many genes that signify pre- 
B cell differentiation, including key genes implicated in IgK gene 
recombination. 

Progressively Diminished and GLTs in Btk~'~ , 
S\p6S~'~, and Btk~'~S\p65~'~ Pre-B Cells 

In these expression profiling studies, we only detected limited 
difierences in germline transcription (GET) over unrearranged Jk 
and Vk gene segments, which is thought to reflect locus 
accessibility [12]. However, we previously showed by serial- 
dilution RT-PCR that the levels of k° 0.8 and k° 1.1 germline 
transcripts, which are initiated in diflerent regions 5'of Jk and 
spliced to the Ck region [49], are apparently normal in Btk ' 
pre-B cells, modestiy reduced in Slp65 ' pre-B cells, and severely 
reduced in Btk~'~ Slp65~'~ pre-B ceUs [34]. We could confirm 
these findings for k" GET by quantitative RT-PCR assays on 
FACS-purified B220"'CD19"' pro-B/pre-B cell fractions 
(Figure 2A). In agreement with our reported findings [34], we 
also found that Btk~^~a.nA Slp65~'~ pre-B cells have defective 
transcription, which is initiated 5' of thejx segments (Figure 2B) 
[49]. 

GET across the Vk region showed a similar pattern of sensitivity 
to pre-BCR signaling: decreased transcription of six individual Vk 
regions tested (Vk3-7, Vk8-24, Vk4-55, Vk10-96, Vk1-35, and 
Vk2-137) correlated with decreased pre-BCR signaling activity 
(Figure 2C) in the pre-B cells of the four groups of mice. GET over 
unrearranged V^l and V\2 segments was strongly reduced in the 
absence of Btk or Slp65, as detected by the expression arrays 
(Table 1). 

These observations indicate that IgK locus accessibilit)', a 
hallmark of recombination-competent antigen receptor loci, is 
progressively reduced under conditions of diminishing pre-BCR 
signaling. 

Pre-BCR Signaling Induces Modulation of Long-Range 
Chromatin Interactions at the Iqk Locus 

Accessibility of antigen receptor loci for V(D)J recombination is 
thought to be initiated by enhancers, in part through long-range 
chromatin interactions with promoters of noncoding transcription, 
resulting in the activation of germline transcription [8] . Because pre- 
BCR signaling affects the expression of GET and various nuclear 
proteins that mediate long-range chromatin interactions and bind 
the K enhancers, it is conceivable that pre-BCR signaling induces 
changes in the enhancer-mediated higher order chromatin structure 
of the IgK locus that facilitates Vk gene accessibility. 

We therefore performed 3C-Seq analyses on FACS-purified 
B220"'"CD19"'" fractions from the same five groups of mice (WT, 
Btk~'~, Slp65~'~, and Btk~'~ Slp65~'~ VH81x transgenic Ragl''' 
pre-B cells, as well as Ragl~'~ pro-B cells). Erythroid progenitors 
were analyzed in pariillcl as a nonlymphoid control, in which the IgK 
locus was not contracted. Genome-wide chromatin interactions were 
measured for three regulatory elements involved in the control of 
locus accessibility and recombination: the IEk and 3'Ek enhancers 
[50-52] and the Sis element [20], which contain binding sites for 
Ikaros/Aiolos, E2a, and Irf4 [16,17,20,38,53]. 

In WT pre-B cells, all three regulatory elements showed 
extensive long-range chromatin interactions within the Vk region 
and substantially less interactions with regions up- or downstream 
of the ~3.2 Mb IgK domain (Figure 3A; see Figure S2, Figure S3, 
and Figure S4 for line graphs), confirming previous observations 
[2 1] . Under conditions of reduced pre-BCR signaling activity, the 
three IgK regulatory elements stiU showed strong interactions with 



PLOS Biology j www.plosbiology.org 



4 



February 2014 | Volume 12 | Issue 2 | e1001791 



Pre-BCR Signaling and IgK Locus Topology 



o 5 



01 _ 

S Q 



■Da 

o O 
u: ^ 



(N ^ ^ 



^ ^ ^ 



00 CX) VO 



^ ^ 



O fN 00 ro a\ 
[N rN LO 0^ 

vd rn rn ro (N 



00 ^ ^ 



^oioq(N rNr^oo^_ 
I I I I ^ I I I I ^ 



in Ln tN (N 



^ Oi CTi VO 



o^^m^ooooa^ m 

I I I I I I I I I ^' ^ ^ I 



ro in ro ro 



m vo 



S o 



iJ2 liH 



< 



< 

Z 
Q 



E E 



< 

Q 



< 
Q 



E E 



o B B 



< 
■z. 

D 



-g s 

Q > 



vD CO 



<Ti ro Ul CO 
"^JOOOOroojcocOvD 



00 vO >- 



i/1 00 pi m 



CO 00 \o 






■T3 ^ 



vD rsi <— ^ 



CM >— 00 T— 



ro >£> Oi vD 



(VI <7i vo ro 
CO \0 O 



\0 ro ro CM 



LnL/imminmi-nmmro 



T- 00 CO 



"= ^ o 
< (ij ^ 



TO 



TO 



> o o 
Q. 



PLOS Biology | www.plosbiology.org 



5 



February 2014 | Volume 12 | Issue 2 | e1001791 



Pre-BCR Signaling and IgK Locus Topology 




K°0.8 



K° 1.1 



B 



<D O 
> in 

,v °- 



I 



A1 




A3 



> m 



V 2-137 

LI 



> m 




V 1-135 




V 4-55 



V 8-24 




Pre-BCR 
signalling 
gradient 




■ WT pre-B 

■ Btk '- pre-B 

□ Slp65-'- pre-B 

□ Btk '-; Slp65-'- pre-B 



Figure 2. Reduction of Btl</Slp65-mediated pre-BCR signaling is associated witii progressive loss of IgK GLT. Quantitative RT-PCR 
analysis for k" (A), X° (B), and GLT (C) of FACS-sorted B220*CD19* pre-B/pro-B cell fractions from the indicated mice on a VH81x transgenic RagV'^ 
background. Expression levels were normalized to those of Gapdh, w/hereby the values in WT pre-B cells were set to one. Bars represent 
mean values and error bars denote standard deviations for four independent mice per group. 
doi:10.1371/journal.pbio.1001791.g002 



the Vk region. Surprisingly, even in tlie complete absence of pre- 
BCR signaling in Ragl~^~ pro-B cells, long-range interactions 
were still observed at frequencies well above those seen in 
nonlymphoid cells, suggesting that a contracted IgK locus 
topology is not strictly dependent on pre-BCR signaling 
(Figure 3A, Figure S2, Figure S3, and Figure S4). Next, we used 
3D DNA FISH analyses using BAC probes hybridizing to the 
distal Vk and CK/enhancer regions to confirm that IgK locus 
contraction was similar in Rag-1 ^ pro-B cells and VH81x 
transgenic Rag-1 ' pre-B cells (both showing a contracted 
topology, compared with noncontracted pre-pro-B cells deficient 
for the TF E2a; Figure 3B). 

Nevertheless, we did observe that pre-BCR signaling induced 
clear diflFerences in interaction frequencies. Whereas an increase in 
pre-BCR signaling was associated with a decrease in the 
interaction frequencies between the two K enhancers and regions 
flanking the IgK locus (as also revealed by more detailed images of 
selected regions upstream and downstream of the IgK domain; see 
Figure S5), the overall interaction frequency within the IgK domain 
appeared unchanged (Figure S3, Figure S4, and Figure S5). 
Remarkably, interactions with the Sis element showed quite an 
opposite pattern: pre-BCR signaling correlated with increased 
overall interactions within the IgK domain and did not substan- 
tially affect interaction frequencies in the IgK flanking regions 
(Figure S2 and Figure S5). 

Taken together, these analyses show that (i) the IgK locus is 
already contracted at the pro-B cell stage and that (ii) pre-BCR 
signaling induces changes in long-range chromatin interactions, 
both within the IgK locus and in the flanking regions. 



Pre-BCR Signaling Enhances Interactions of 3'Ek and Sis, 
But Not IEk, with Fragments 

The differential effects of pre-BCR signaling on long-range 
chromatin interactions of the IEk, 3'Ek, and Sis elements clearly 
emerged in a quantitative analysis of the 3C-seq datasets (Figure 4A; 
see Materials and Methods for a detailed description of the 
quantification methods used). When pre-BCR signaling was absent 
{Ragl~'~ pro-B cells) or very low [BtUT ' ~ Slp65~ ' ~ pre-B cells), the 
average interaction frequencies were similar within the ~3.2 Mb 
Vk region and the ~3.2 Mb downstream flanking region, for all 
three regulatory elements. Interaction frequencies with the 
upstream flanking region were lower, consistent with the larger 
chromosomal distance to the three viewpoints. The presence of 
increasing levels of Btk/Slp65-mediated pre-BCR signaling was 
associated with reduced interaction of IEk and 3'Ek with the IgK 
flanking regions and with increased interaction of the Sis element 
and (to a lesser extent) 3'Ek with the Vk region (Figure 4A). As a 
result, for all three regulatory elements pre-BCR signaling resulted 
in a preference for interaction with fragments inside the Vk region 
over fragments outside the Vk region (Figure S7). 

We next focused our analysis on the Vk region and compared 
fragments that harbor a functional Vk gene (Vk^ fragment) and those 
that do not (Vk~ fragment). When pre-BCR signaling was absent 
[Ragl ' pro-B cells) or very low iji^ ' Slp65 ^ pre-B cells), the 
average interaction frequencies of the Sis or IEk elements with Vk^ 
fragments were higher than with Vk~ fragments. The average 
interaction frequencies of 3'Ek with Vk^ and Vk~ fragments, 
however, were similar (Figure 4B). Upon pre-BCR signaling, the Sis 
element showed an increase in interaction frequencies with both Vk^ 
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Figure 3. 3C-Seq analysis of iong-range cliromatin interactions witKiin tiie IgK locus and flanking regions. (A) Overview of long-range 
interactions revealed by 3C-Seq experiments performed on the indicated cell fractions, representing a gradient of pre-BCR signaling, whereby the IEk 
element (top), the 3'Ek element (center), or the Sis element (bottom) was used as a viewpoint. Shown are the relative interaction frequencies 
(average of two replicate experiments) for the indicated genomic locations. The —8.4 Mb region containing the IgK locus (yellow shading) and 
flanking regions (cyan shading) is shown and genes and genomic coordinates are given (bottom). The locations of the two BAC probes used for 3D 
DNA-FISH are indicated by a green (distal Vk probe) and red (proximal Cx/enhancer probe) rectangle {bottom). Pre-B cell fractions were FACS-purifled 
from the indicated mice on a VH81x transgenic Ragl~'~ background (see Figure 1 for gating strategy). Erythroid progenitor cells were used as a 
nonlymphoid control. (B) 3D DNA-FISH analysis comparing locus contraction in cultured bone-marrow-derived E2a~'~ pre-pro-B, Ragl~'~ pro-B, and 
VI-I81X Ragl~'~ pre-B cells (see Figure S6 for phenotype of IL-7 cultured B-lineage cells). Locations of the BAC probes used are indicated at the 
bottom of panel A. Representative images for all three cell types are shown on the left, quantifications (>100 nuclei counted per cell type) on the 
right. The red lines indicate the median distance between the two probes. Statistical significance was determined using a Mann-Whitney U test 
(***p<0.001; n.s., not significant, paO.OS). 
doi:l 0.1 371/journal.pbro.lOOl 791 .g003 



and Vk~ fragments, with nevertheless an interaction preference 
for Vk^ fragments. In contrast, interaction frequencies between 
the iEk element and or Vk~ fragments were not modulated 
by pre-BCR signahng at all (Figure 4B). The 3'Ek element 
exhibited yet another profile: pre-BCR signaling induced 
increased interaction frequencies specifically with Vk^ fragments, 
while interactions with Vk~ fragments were not notably 
modulated by pre-BCR signaling (Figure 4B). When we 
separately analyzed nonfunctional pseudo-VK genes, we found 
for the Sis and 3'Ek elements that the interaction patterns with 
functional and nonfunctional Vk genes were similar (Figure S8). 
In contrast, the IEk enhancer did show an overall increased 
interaction frequency with Vk functional genes, compared with 
nonfunctional Vk genes, a phenomenon which was again 
independent from pre-BCR signaling (Figure S8). 

The finding that interactions of Vk genes with the intronic 
enhancer are already robust in pro-B cells, while those with the 
3'k enhancer are dependent on pre-BCR signaling, suggested that 
for individual Vk genes pre-BCR signaling may result in more 
similar interaction frequencies with the two enhancers. To 
investigate this, we examined for all individual Vk genes the 
correlation between their 3C-seq interaction frequencies with the 
IEk and 3'k elements and found that these were highly correlated 
in W'T pre-B cells (R^ = 0.68; Figure 4C). Correlation was severely 
reduced when pre-BCR signaling was low in Btk ' Slp65 ' pre- 
B ceUs (R^ = 0.26; Figure 4C). Similar pre-BCR signaling- 
dependent correlations were observed between VK-interactions 
with the Sis element and those with the two enhancers (Figure S9). 
As the Sis element particularly suppresses recombination of the 
proximal Vk3 family, we investigated interaction correlations 
specifically for this Vk family. Similar to our findings for aU Vk 
genes, a subanalysis showed strong correlations for the interactions 
of Vk3 family genes with IEk, 3'k, and Sis in WT pre-B cells, 
which \v(;re diminished when prc'-BCR signaling w^Vn lo\v, except 
for iEK-Sis correlations, which were pre-BCR signaling-indepen- 
dent (Figure S9). 

In summary, we conclude that pre-BCR signaling induces a 
redistribution of long-range interactions of the iEK, 3'Ek, and 
Sis elements, tlicrcljy restricting interactions towards the V^ 
gene region. Moreover, upon pre-BCR signaling the long- 
range interactions mediated by 3'Ek and Sis — but not those 
mediated by iEk — become enriched for fragments harboring a 
Vk gene, demonstrating increased proximity of 3'Ek and Sis 
to Vk genes. Finally, for individual Vk genes, the interactions 
with IEk, 3'Ek, and Sis become highly correlated upon pre- 
BCR signaling, indicating that pre-BCR signals result in 
regulatory coordination between these three elements that 
govern IgK locus recombination. In contrast, interactions 
between genes of the proximal Vk3 family. Sis and IEk — but 
not 3'k — appear to be coordinated already in the absence of 
pre-BCR signaling. 



Long-Range Chromatin Interactions of k Regulatory 
Elements Correlate with Vk Usage 

Next, we investigated the effects of pre-BCR signaling on the 
interaction frequencies of individual functional Vk genes with the 
three K regulatory elements (Figure 5A,B). The 3C-seq patterns of 
the majority (~91%) of the 101 individual Vk^ fragments showed 
evidence for interaction with one or more of the K regulatory 
elements (> 25 average counts). When comparing _SiA ' Slp65 ' 
with WT pre-B cells, we observed that for a large proportion 
(~38-52'^i)) of Vk fragments, interaction frequencies increased 
upon pre-BCR signaling (Figure 5B). Smaller proportions of Vk^ 
fragments showed a decrease (~ 12-29%) or were not significandy 
affected by pre-BCR signaling (~ 17-25% with < 1.5-fold change). 
The observed increase or decrease was not related to proximal or 
distal location of the Vk genes, nor to their sense or antisense 
orientation (not shown). Distributions of the three different classes 
of Vk fragments showed substantial differences between the K 
regulatory elements. For the Sis and 3'Ek elements, more Vk^ 
fragments showed increased than decreased interactions 
(Figure 5B), in agreement with the signaling-dependent increase 
in average interaction frequencies of all Vk^ fragments (Figure 4B). 
In contrast, for the IEk viewpoint, Vk^ fragments showing 
increased and decreased interactions were more equal in number, 
consistent with the limited effects of pre-BCR signaling on overall 
iEK interaction frequencies of aU Vk* fragments (Figure 4B). 

Although antigen receptor recombination is in principle regarded 
as a random process, a significant skewing of the primary IgK 
repertoire of C57BL/ 6 mice was recendy reported: one third of the 
Vk genes was shown to account for >85"/o of the Vk segments used 
by B cells [54] . To assess whether a correlation exists between usage 
of Vk genes and their interaction frequencies with K regulatory 
elements, we divided the Vk genes into four usage categories (< 
0.1%, 0.1-0.3%, 0.3-0.5%, and >0.5%) and calculated their 
average 3C-Seq interaction frequencies with Sis, IEk, and 3'k 
(Figure 5C). In WT pre-B cells, Vk usage showed a strong positive 
correlation with 3C-Seq interaction frequencies for all three 
regulatory elements (R^ = ~0.7-0.9; Figure 5C). These correlations 
were pre-BCR signaling-dependent, since in Bt^' Slp65~' pre-B 
cells, they were reduced (for IEk; R^ = 0.33) or absent (for Sis and 
3'k; r2<0.10 and R^<0.16, respectively) (Figure 5C). 

Collectively, our results indicate that .specifically the most frequently 
used Vk genes are the main interaction targets of K regulatory 
elements, whereby pre-BCR signaling completely underlies this 
specificity for the Sis and 3'Ek elements, and to a lesser extent for iEk. 

Long-Range Interactions with k Regulatory Elements 
Correlate with the Presence of Ctcf, Ikaros, E2a, and H3K4 
Hypermethylation 

Next, we investigated whether long-range interactions between 
K regulatory elements and the Vk region correlated with the 
presence of the TFs Ctcf [21], Ikaros [55], and E2a [56], which 
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Figure 4. Modulation of long-range chromatin interactions within the IgK locus by pre-BCR signaling. Quantitative analysis of 3C-Seq 
datasets using the three indicated k regulatory elements as viewpoints. (A) Average long-range chromatin interaction frequencies (from two replicate 
3C-seq experiments) with upstream (—2.0 Mb), Vk (—3.2 Mb), and downstream (—3.2 Mb) regions, as defined in Figure 3A, for the five B-cell 
precursor fractions representing a pre-BCR signaling gradient. Average interaction frequencies per region were calculated as the average number of 
3C-Seq reads per restriction fragment within that region. See Materials and Methods section for more details. (B) Average interaction frequencies 
within the Vk region were determined for fragments that do not (— ) contain a functional gene and for those that do contain a functional Vk gene 
(+). (C) Correlation plots of average interaction frequencies of the two enhancer elements with the 101 functional V^ genes for WT pre-B cells (left) 
versus Btk~'~Slp65~'~ pre-B cells (right). On the log scale, frequencies <1 were set to 10°. Statistical significance was determined using a Mann- 
Whitney U test (*p<0.05; **p<0.01; ***p<0.001; n.s., not significant, p>0.05). 
doi:1 0.1 371/journal.pbio.1 001 791 .g004 



have been implicated in IgK locus recombination [21,37,55,57,58]. 
Notably, Ikaros and E2a both strongly bind all three K regulatory 
elements, while the Sis element is also occupied by Ctcf ([21]; 
unpublished data). 

Remarkably, we found similar striking correlations between the 
presence of in vivo binding sites for each of these TFs (as 
determined by ChIP experiments; see Materials and Methods for 
the relevant references) and long-range chromatin interactions 
with the K regulatory elements (Figure 6A-C), even though Ctcf 
sites are mostly located in between Vk genes [2 1] and Ikaros/E2a 
sites were frequently found close to Vk gene promoter regions ([2]; 
Figure 7A). Even when pre-BCR signaling was absent (Ragl~'~ 
pro B cells) or very low (Btk ' Slp65 ' pre-B cells), the average 
interaction frequencies of the K regulatory elements with fragments 
containing Ctcf, Ikaros, or E2a bindings sites were higher than 
those without binding sites. Irrespective of the presence or absence 
of bindings sites for these TFs, we found that upon pre-BCR 
signahng interaction frequencies with the Sis element increased 
and those with the iEk did not change. In contrast, for the 3'Ek 
we found that pre-BCR signaling specifically increased interaction 
frequencies with fragments occupied by Ctcf, Ikaros, or E2a. 

Finally, we found that the presence of di- or trimethylation of 
histone 3 lysine 4 (H3K4Me2/3), an epigenetic signature associated 
with locus accessibilit}' [59] and Rag-binding [60,61], also 
correlated with increased interaction frequencies with K regulatory 
elements, revealing a similar pre-BCR signaling dependency as seen 
for the TFs analyzed (Figure 6D). 

We conclude that the presence of essential TFs or H3K4Me2/3 
in the Vk region strongly correlates with the formation of long- 
range chromatin interactions with the K regulatory elements, and 
that for the Sis and 3'Ek elements this interaction preference is 
further enhanced by pre-BCR signaling. 

Proximity of Vk Genes to E2a Binding Sites Correlates 
with Higli Vk Usage and Increased Long-Range 
Chromatin Interactions 

Since the long-range interactions with K regulatory elements 
correlated with the presence of TFs implicated in IgK. recombina- 
tion, we next asked whether the k regulatory elements preferen- 
tially interacted with Vk genes that are in close proximity to 
binding sites for Ctcf, Ikaros, or E2a. 

Strikingly, the majority of functional Vk genes (95/101) was 
found to have an Ikaros binding site in close proximity — that is, 
located on the same 3C-seq restriction fragment (average length of 
~3 kb, unpublished data) (Figure 7A). Proximity of Vk genes to an 
E2a binding site (37%) or H3K4Me2/3 positive region (~28%) is 
more selective, while only a small fraction of Vk genes are c:lose to 
Ctcf binding sites (~ 1 2"/a) ([22]; Figure 7A). All Vk genes marked by 
E2a, Ctcf, H3K4Me2/3, or a combination of these also contain an 
Ikaros binding site. Frequently used Vk genes (>1.0% usage; 33/ 
101 genes) were located in two separate regions, a proximal and a 
distal region, which also contained virtually all E2a and H2K4Me2/ 
3-marked Vk genes (Figure 7A). 



We found that Vk genes marked by both Ikaros and E2a were 
used substantially more often than those only bound by Ikaros 
(Figure 7B), suggesting that these Vk genes are preferentially 
targeted for Vk-Io-Jk gene rearrangement. Our 3C-seq analyses 
showed that in WT pre-B cells, interaction frequencies with the 
three k regulatory elements were higher for Ikaros/E2a-marked 
Vk genes compared to genes marked by Ikaros binding alone 
(Figure 7C). In fact, Vk^ restriction fragments containing an Ikaros 
binding site but not an E2a binding site showed interaction 
frequencies similar to Vk~ restriction fragments. Under conditions 
of very low pre-BCR signaling (in Btk ' Slp65 ' pre-B cells), we 
observed strongly reduced interaction frequencies of Vk^ E2a 
binding restriction fragments with the Sis and 3'Ek elements. 
These interaction frequencies were in the same range as those of 
Vk fragments or Vk fragments that harbored an Ikaros site only 
(Figure 7C). Interaction frequencies with the IEk enhancer, 
however, were independent of pre-BCR signaling. As shown in 
Figure 7D, for the majority of Ikaros/E2a-marked Vk^ fragments 
(65%), pre-BCR signaling was associated with increased interac- 
tions with the Sis and 3'Ek elements (comparing wild-type and 
Blk~^~ Slp65~^~ pre-B cells). In these analyses, only ~13.5% and 
~5.4% of Ikaros/E2a-marked Vk^ fragments showed a decreased 
interaction frequency upon pre-BCR signaling. In contrast, almost 
equal proportions of Ikaros/E2a-marked Vk^ fragments showed 
increased (~37%) and decreased (~30%) interactions with iEK 
upon pre-BCR signaling. 

Taken together, these data reveal strong positi\'e correlations 
between the presence of E2a binding sites, Vk usage, and long- 
range chromatin interactions with K regulatory elements in pre-B 
cells. Remarkably, for the IEk element, these correlations are 
largely independent of Btk/Slp65-mediated pre-BCR signaling, 
whereas for the 3'Ek they are completely dependent on signaling. 

Discussion 

During B-cell development the pre-BCR checkpoint is known to 
regulate the expression of many genes, part of which control the 
increase in IgK locus accessibility to the V(D}J recombinase 
complex. However, it remained unknown how pre-BCR signaling 
events affect accessibility in terms of IgK locus contraction and 
topology. 

Here we identified numerous genes involved in IgL chain 
recombination, diromatin modification, signaling, and cell 
survival to be aberrandy expressed in pre-B cells lacking the pre- 
BCR signaling molecules Btk and/or Slp65. We found that GLT 
over the Vk region, reflecting Vk accessibility, is strongly reduced 
in these cells. We used 3C-Seq to show that in pro-B cells both the 
intronic and the 3' K enhancers frequendy interact with the 
~3.2 Mb Vk region, as well as with IgK flanking sequences, 
indicating that the IgK locus is already contracted at the pro-B cell 
stage. 3C-Seq analyses in wild-t)'pe and Btk/Slp65 single- and 
double-deficient pre-B cells demonstrated that pre-BCR signahng 
significandy affects IgK locus topology. First, pre-BCR signaling 
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Figure 5. Long-range chromatin interactions of k regulatory elements correlate with Vk gene usage. (A) Selected examples of genomic 
regions containing V/ fragments, showing increased (Vk4-57, \/„6-25, \/„8-24, \/„4-79), stable (\/„8-23-^, \/^4-78, Vk13-84), or decreased (V^S-Q, 
V^14-111) 3C-seq interaction frequencies with 3'Ek or iEic upon pre-BCR signaling. Averaged 3C-seq signals are plotted as a line graph, with the 
individual data points representing the center of the Bglll restriction fragments. Yellow shading marks the Bglll fragment on which the gene(s) is 
located. V,^ gene(s) are indicated (top) and chromosomal coordinates and scale bars (10 kb) are plotted (bottom). (B) Classification of V^* fragments, 
based on the effect of pre-BCR signaling on their interactions with the three k: regulatory elements indicated. Increase and decrease were defined as 
>1. 5-fold change of interaction frequencies detected in WT pre-B cells versus Btk^'^SIp65^'^ pre-B cells. (C) Correlation of average interaction 
frequencies (for the three k: regulatory elements indicated) with four usage categories ranging from low (<0.1%) to high usage (>0.5%, listed in 
the table on the right). Diamonds represent average interaction frequencies for Btk^'^SIp65^'^ pre-B cells (yellow) and WT pre-B cells (grey). The 
dotted line in each graph depicts the average interaction frequency with fragments that do not contain a functional (V,^-^). Primary V,^- gene usage 
data were taken from [54]. 
doi:10.1371/journal.pbio.1001791.g005 



reduces the interactions of the intronic and 3 ' K enhancers with IgK 
flanking regions, efFectively focusing enhancer action towards the 
Vk region to facilitate Vk-Io-Jk recombination. Second, pre-BCR 
signaling strongly increases nuclear proximity of the 3'k enhancer 



to Vk genes, whereby this increase is more substantial for more 
frequently used Vk genes and for Vk genes close to a binding site 
for the basic helix-loop-helix protein E2a. Third, pre-BCR 
signaling augments interactions between K regulatory elements 
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Figure 6. Long-range chromatin interactions of k regulatory elements correlate with TF binding and histone modifications. (A-D) For 
fragments within the region, average 3C-seq interaction frequencies were calculated for fragments that did (+) or did not (-) contain binding sites 
for TFs or H3K4 histone modifications (as determined by previous ChlP-Seq studies; see IVIaterials and Methods for references). Data for the three 
viewpoint and the five B-cell precursor fractions representing a pre-BCR signaling gradient are shown for Ctcf (A), Ikaros (B), E2a (C), and H3K4 di- and 
tri-methylation (Me2/3). Statistical significance was determined using a IVIann-Whitney U test (*p<0.05; **p<0.01; ***p<0.001; n.s., not significant, 
p>0.05). 

doi:l 0.1 371/journal.pbro.lOOl 791 .g006 



and fragments within the region bound by the key B-cell TFs 
Ikaros and E2a and the architectural protein Ctcf. Fourth, pre- 
BCR signaling has limited efTects on interactions of the intronic K 
enhancer with fragments within the IgK locus, as this enhancer 
already displays interaction specificity for functional Vk genes and 
TF-bound regions in pro-B cells. Fifth, pre-BCR signaling has 
limited effects on the interactions between the intronic or 3'k 
enhancers and fragments that do not contain a gene or an 
Ikaros, E2a, or Ctcf binding site, emphasizing the specificity of 
pre-BCR signaling-induced changes in IgK locus topology. Sbcth, 
pre-BCR signaling appears to induce mutual regulatory coordi- 
nation between the three regulatory elements, as their interaction 
profiles with individual Vk genes become highly correlated upon 
signaling. Finally, pre-BCR signaling increases interactions of the 
Sis element with DNA fragments in the IgK locus, irrespective of 
the presence of a Vk gene or TF. Collectively, our findings 
demonstrate that pre-BCR signals relayed through Btk and SlpG.T 
are required to create a chromatin environment that facilitates 
proper IgK locus recombination. This multistep process is initiated 
by up-regulation of key TFs like Aiolos, Ikaros, Irf4, and E2a. 
These proteins are then recruited to or further accumulate at the 
IgK locus and its regulatory elements, resulting in a specific fine- 
tuning of enhancer-mediated locus topology that increases locus 
accessibility to the Rag recombinase proteins. 

Importandy, the presence of strong lineage-specific interaction 
signals between the CK/enhancer region and distal Vk genes in 
pro-B cells indicates that the IgK locus is already contracted at this 
stage. In contrast to a previous microscopy study indicating that 
IgK locus contraction did not occur until the small pre-B cell stage 
[36], our 3D DNA FISH analysis indeed detected similar nuclear 
distances between distal Vk and the CK/enhancer region in 
cultured pro-B and pre-B cells. Recently Hi-C was employed to 
study global early B cell genomic organization whereby substantial 
interaction frequencies were found between the intronic K 
enhancer and the Vk region in pro-B cells [40]. E2a-deficient 
pre-pro-B cells, which are not yet fuUy committed to the B-ceU 
fineage [62], showed very few interactions among the IEk and the 
distal part of the Vk region [40], resembling the interactions we 
observed in nonlymphoid cells (Figure 3A). Accordingly, 3D-FISH 
analysis showed that the IgK locus adopted a noncontracted 
topology in these pre-pro-B cells (Figure 3B). These data indicate 
that IgK locus contraction is already achieved in pro-B cells and 
depends on the presence of E2a. Supporting this notion, active 
histone modifications and E2a were already detected at the K 
enhancers and Vk genes at the pro-B cell stage [.56,63], whereby 
E2a was frequendy found at the base of long-range chromatin 
interactions together with Ctcf and Pu. 1, possibly acting as 
"anchors" to organize genome topology [40]. The observed 
correlation between E2a binding, Vk gene usage and iEk 
proximity in pro-B cells (Figure 5C, Figure 7C) further strengthens 
an early critical role for E2a in regulating IgK locus topology, Vk 
gene accessibility, and recombination. 

Our 3C-seq experiments revealed that pre-BCR signaling is not 
required to induce long-range interactions between the K 
regulatory elements and distal parts of the Vk locus, indicating 
that TFs strongly induced by signaling — that is, Aiolos, Ikaros, and 



Irf4 — are not strictly necessary to form a contracted IgK locus. 
Prime candidates for achieving IgK locus contraction at tlu- pro-B 
cell stage are E2a and Ctcf, as they have been implicated in 
regulating Ig locus topology [21,40,64,65] and E2a already marks 
frequently used Vk genes at the pro-B cell stage (Figure 7), 
although we did obserx^e reduced E2a expression and binding to 
the iEK enhancer and Vk genes when pre-B cell signaling was low 
(Figure 1 and Table S3), suggesting that pre-BCR signaling is 
required for high-level E2a occupancy of the Vk genes. We 
previously reported that IgK gene recombination can occur in the 
absence of Ctcf and that Ctcf mainly functions to limit interactions 
of the K enhancers with proximal Vk regions and to prevent 
inappropriate interactions between these strong enhancers and 
elements outside the IgK locus [21]. Because at the pro-to-pre-B 
cell transition Aiolos, Ikaros, and Irf4 are recruited to the IgK locus 
and histone acetylation and H3K4 methylation increases 
[17,38,63,66], we hypothesize that pre-BCR-induced TFs act 
upon an E2a/Ctcf-mediated topological scaffold to further refine 
the long-range chromatin interactions of the k regulatory 
elements. Hereby, these TFs mainly act to focus and to coordinate 
the interactions of the two K enhancers to the Vk gene s(;gmc'nts, in 
particular to frequently used Vk genes, thereby increasing their 
a[:cessibility for recombination (see Figure 7E for a model of pre- 
BCR signaling-induced changes in IgK locus accessibility). 

In this context, our 3C-seq data show that the two k enhancer 
elements have distinct roles. Both 3'Ek and IEk elements manifest 
interaction specificity for highly used, E2a-marked, Vk genes. 
However, whereas IEk already shows this specificity in pro-B cells 
(although pre-BCR signaling does augment this specificity), 3'Ek 
only does so in pre-B cells upon pre-BCR signaling. These 
observations indicate that IEk is already "prefocused" at the pro-B 
cell stage and that pre-BCR signals are required to fully activate 
and focus the 3'Ek to allow synergistic promotion of IgK 
recombination by both enhancers (see Figure 7E) [52]. In 
agreement with such distinct sequential roles, IEk and not the 
3'Ek was found to be required for the initial increase in IgK locus 
accessibility, which occurred upon binding of E2a only [37,38,67]. 
The 3'Ek on the other hand requires binding of pre-BCR 
signaling-induced Irf4 to promote locus accessibility [19,38], 
followed by further recruitment of E2a to both k enhancers and 
highly used Vk genes (Table S3 and [38,57]). 

The Sis regulatory element was shown to dampen proximal Vk— 
Jk rearrangements and to specify the targeting of .^K transgenes to 
centromeric heterochromatin in pre-B cells [20]. As Sis is 
extensively occupied by the architectural Ctcf protein and deletion 
of Sis or Ctcf both resulted in increased proximal Vk usage 
[21,23], it was postulated that Sis functions as a barrier element to 
prevent the k enhancers from too frecjuently targeting proximal Vk 
genes for recombination. In this context, we now provide evidence 
that interactions between the proximal Vk genes. Sis, and IEk — 
but not 3'k — are already coordinated before pre-BCR signaling 
occurs (Figure S9). Perhaps not surprisingly. Sis-mediated long- 
range chromatin interactions displayed a pattern and pre-BCR 
signaling response that was different from the K enhancers. Unlike 
for the enhancers, upon pre-BCR signaling. Sis-mediated interac- 
tions with regions outside the IgK locus were maintained and 
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Figure 7. Proximity of genes to E2a binding sites correlates witli frequencies of long-range interactions. (A) Schematic 
representation of the Iqk locus, showing the location of all functional (grey, top), J.^ and Q gene segments, and the k regulatory elements Sis, IEk, 
and 3'Ek. IVIAR, matrix attachment region. Vk genes within close proximity (as defined by colocalization on the same 3C-Seq restriction fragment) to 
the indicated TPs or H3K4 hypermethylation (as detected by previous ChlP-seq studies; see IVlaterials and Methods for references) are shown. At the 
bottom, highly used (>1.0% used) Vk gene segments are depicted (orange), which cluster within two large high-usage domains (yellow shading). 
Primary Vk gene usage data was taken from [54]. (B) Average usage of Vk genes marked only by an Ikaros binding site or those marked by binding 
sites of both Ikaros and E2a. (C) Comparison of average interaction frequencies (for the three k regulatory elements indicated) between Vk~ 
fragments (no VJ, Vk* fragments containing an Ikaros binding site only, and Vk* fragments containing both an Ikaros and E2a binding site. Bars 
represent average frequencies for Btk~'~Slp65~'~ pre-B cells (yellow) and WT pre-B cells (grey). (D) Classification of Vk* fragments, containing an 
Ikaros binding site only (top) or containing both an Ikaros and E2a binding site (bottom), based on the effect of pre-BCR signaling on their 
interactions with the three k regulatory elements indicated. Increase and decrease were defined as >1. 5-fold change of interaction frequencies 
detected in WT pre-B cells versus Btk~'~Slp65~'~ pre-B cells. (E) Proposed model of pre-BCR signaling-mediated changes in k enhancer action. In pro- 
B cells (left) the enhancers show minimal coordination and their interactions are not yet (fully) focused on the Vk genes. Upon pre-BCR signaling and 
differentiation to pre-B cells (right), TPs bind the locus to coordinate enhancer action and focus their interactions to the Vk genes, inducing germline 
transcription (GLT) and accessibility to the V(D)J recombinase. See Discussion for more details. Statistical significance was determined using a IMann- 
Whitney U test (*p<0.05; **p<0.01; ***p<0.001; n.s., not significant, p>0.05). 
dor:l 0.1 371/journal.pbio.lOOl 791 .g007 



interaction within the Vk region increased, irrespective of the 
presence of Vk genes or TF binding sites. Because Sis is involved in 
targeting the nonrecombining IgK allele to heterochromatin [20], 
the observed interaction pattern of the Sis element might reflect its 
action in pre-B cells to sequester the nonrecombining IgK locus 
and target it towards heterochromatin. This might also explain the 
increased interaction frequencies of Sis with highly used Vk genes 
upon pre-BCR signaling (Figures 5C and 7C), as such highly 
accessible genes likely require an even tighter association with Sis 
and heterochromatin to prevent undue recombination. 

Surprisingly, we observed a striking correlation between Ikaros 
binding and Vk gene location (94% of Vk genes were in close 
proximity to an Ikaros binding site; Figure 7A). Although Ikaros 
and Aiolos have a positive role in regulating gene expression 
during B-cell development [55,58] and Ikaros is required for IgH 
and IgL recombination [39,58], Ikaros has also been reported to 
silence gene expression through its association with pericentro- 
meric heterochromatin [68] or through recruitment of repressive 
cofactor complexes [69,70]. R(;(Tuitmc-nt of Ikaros to the IgK locus 
was found increased in pre-B cells as compared to pro-B cells [63], 
in agreement with its up-regulation in pre-B cells (Figure 1). 
Furthermore, Ikaros binds the Sis clement, where it was suggested 
to mediate heterochromatin targeting of IgK alleles by the Sis 
region [20]. Aiolos, although not essential for B-cell development 
like Ikaros [58,71], is strongly induced by pre-B cell signaling and 
has been reported to cooperate with Ikaros in regulation gene 
expression [27]. Although their synergistic role during IgL chain 
recombination has not been extensively studied, the Ikaros/Aiolos 
ratio changes upon pre-BCR signaling (Figure 1). Increased 
recruitment of Ikaros/Aiolos to Vk genes and the K enhancers 
likely increases IgK locus accessibility and contraction (see 
Figure 6), as Ikaros was very recendy shown to be essential for 
IgL recombination [58] . On the other hand, it is conceivable that 
on the nonrecombining allele, increased recruitment of Ikaros/ 
Aiolos to Vk genes and the Sis region could facilitate silencing of 
this allele. Further investigations using alkJe-specific approaches 
[72] win be required to clarify the allele-specific action of the Sis 
element during IgK recombination. 

In summary, by investigating the effects of a pre-BCR signaling 
gradient — rather than deleting individual TFs — ^we have taken a 
more integrative approach to study the regulation of IgK locus 
topology. Our 3C-Seq analyses in wild-type, Btk, and Slp65 single- 
and double-deficient pre-B cells show that interaction frequencies 
between Sis, IEk, or 3' Ek and the Vk region are already high in 
pro-B cells and that pre-BCR signaling induces accessibility 
through a functional redistribution of long-range chromatin 
interactions within the Vk region, whereby the IEk and 3'Ek 
enhancer elements play distinct roles. 



Materials and Methods 

Mice 

VH81x transgenic mice [73] on the Rag-1~^~ background [74] 
that were either wild-type, Btk~'~ [75], Slp65~'~ [42], or BtT' 
~Slp65~'~ have been previously described [34]. Mice were 
crossed on the C57BL/6 background for >8 generations, bred, 
and maintained in the Erasmus MC animal care facility under 
specific pathogen-free conditions and were used at 6-13 wk of age. 
Experimental procedures were reviewed and approved by the 
Erasmus University Committee of Animal Experiments. 

Flow Cytometry 

Preparation of single-cell suspensions and incubations with 
monoclonal antibodies (mAbs) were performed using standard 
procedures. Bone marrow B-lineage cells were purified using 
fluorescein isothiocyanate (FITC)-conjugated anti-B220(RA3-6B2) 
and peridinin chlorophyll protein (PCP)-conjugated anti-CD 19, 
together with biotinylated mAbs specific for lineage markers Gr- 1 , 
Terl 19, and CDl lb and APC-conjugated streptavidin as a second 
step to further exclude non-B cells. Cells were sorted with a 
FACSARia (BD Biosciences). The following mAbs were used for 
flow cytometry: FITC-, PerCP-anti-B220 (RA3-6B2), phycoery- 
thrin (PE)-anti-CD2 (LFA-2), PCP-, allophycocyanin (APC)- or 
APC-Cy7-anti-CD19 (ID3), PE-, or APC anti-CD43 (S7). All 
these antibodies were purchased from BD Biosciences or 
eBiosciences. Samples were acquired on an LSRII flow cytometer 
(BD Biosciences) and analyzed with Flowjo (Tree Star) and 
FACSDiva (BD Biosciences) software. 

Quantitative RT-PCR and DNA Microarray Analysis 

Extraction of total RNA, reverse-transcription procedures, 
design of primers, and cDNA amplification have been described 
previously [21]. Gene expression was analyzed using an ABI Prism 
7300 Sequence Detector and ABI Prism Sequence Detection 
Software version 1.4 (Apphed Biosystems). AU PCR primers used 
for quantitative RT-PCR of TFs or k°, A,°, and Vk GLT are 
described in [21], except for Obfl (forward 5'-CCTGGCCACC- 
TACAGCAC-3', reverse 5'-GTGGAAGCAGAAA CCTCCAT- 
3', obtained from the Roche Universal Probe Library). 

Biotin-labeled cRNA was hybridized to the Mouse Gene 1 .0 ST 
Array according to the manufacturer's instructions (Affymetrix); 
data were analyzed with BRB-ArrayTools (version 3.7.0, National 
Cancer Institute) using Afiymetrix CEL files obtained from GCOS 
(Affymetrix). The RMA approach was used for normalization. The 
TIGR MultiExperiment Viewer software package (MeV version 
4.8.1) was used to perform data analysis and visuahze results [45]. 
One-way ANOVA analysis of the five experimental groups of B 
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cells was used to identify genes significantly different from wild- 
type VH81X Tg Ragr'~ pre-B ceUs {p<Q.Q\). 

Chromatin Immunoprecipitation (ChIP) 

Chip experiments were performed as previously described [76] 
using FACS sorted bone marrow pre-B cell fractions (0.3-2.0 
miUion cells per ChIP). Antibodies against E2a (sc-349, Santa 
Cruz Biotechnology) and Ikaros (sc-9861, Santa Cruz Biotechnol- 
ogy) were used for immunoprecipitation. Purifu'd DNA was 
analyzed by quantitative RT-PCR as described above. Primer 
sequences are available on request. 

Chromosome Conformation Capture Coupled to High- 
Throughput Sequencing (3C-Seq) 

3C-Seq experiments were essentially carried out as described 
previously [21,41]. For 3C-Seq library preparation, Bglll was 

used as the primary restriction enzyme and Nlalll a.s a secondary 
restriction enzyme. 3C-seq template was prepared from WT 
E13.5 fetal liver erythroid progenitors and FACS-sorted bone 
marrow pro-B cell or pre-B cell fractions (see above) from pools of 
4—6 mice. In total, between 1 and 8 million cells were used for 
3C-seq analysis. Primers for the Sis, IEk, and 3'Ek viewpoint- 
specific inverse PCR were described previously [21]. 3C-seq 
libraries were sequenced on an lUumina Hi-Seq 2000 platform. 
3C-Seq data processing was performed as destxibed elsewhere 
[41,77]. Two repKcate experiments were sequenced for each 
genotype and viewpoint, and normalized interaction frequencies 
per Bgill restriction fragment were averaged between the two 
experiments. 

For quantitative analysis, the IgK locus and surrounding 
sequences were divided into three parts (mm9 genome build): a 
~2 Mb upstream region (chr6:65,441, 978-67,443,029; 759 frag- 
ments), a -3.2 Mb Vk region (chr6:67,443,034-70,80 1,754; 1,290 
fragments) and a downstream ~3.2 Mb region (chr6:70,801,759- 
73,993,074; 1,143 fragments). For each cell type (as described 
above) sequence read counts within individual Bglll restriction 
fragments were normalized for differences in library size 
(expressed as "reads per million"; see [74]) and averaged between 
the two replicates before further use in the various calculations. 
Very small Bglll fragments (<100 bp) were excluded from the 
analysis. Fragments in the immediate vicinity of the regulatory 
elements (chr6:70,659,392-70,693,183; 10 fragments) were also 
excluded because of high levels of noise around the viewpoint, a 
characteristic of all 3C-based experiments. Vk gene coordinates 
(both functional genes and pseudogenes) were obtained from 
IMGT [1 1] and NCBI (Gene ID: 243469) databases. Vk gene 
usage data (C57BL/6 strain, bone marrow) were obtained from 
[54]. ChlP-seq datasets were obtained from [21] (Ctcf), [55] 
(Ikaros), and [56] (E2a, H3K4Me2, and H3K4Me3). Vk genes 
were scored positive for TF binding sites or for a histone 
modification, if they were located on the same Bglll restriction 
fragment (corresponding to the 3C-Seq analysis). 

3D DNA Immuno-FISH 

Kag-l pro- and Rag- 1 ' ;VH8 1 X pre-B ceUs were isolated 
from femoral bone marrow suspensions by positive enrichment of 
CD19^ cells using magnetic separation (MUtenyi Biotec). Cells 
were cultured for 2 wk in Iscove's Modified Dulbecco's medium 
containing 10% fetal calf serum, 200 U/ml penicillin, 200 mg/ml 
streptomycin, 4 nM L-glutamine, and 50 \xM f3-mercaptoethanol, 
supplemented with IL-7 and stem cell factor at 2 ng/ml. E2a~'~ 
hematopoietic progenitors were grown as described previously 
[78] . Prior to 3D-FISH analysis, cells were characterized by flow 



cytometric analysis of CD43, CD 19, and CD2 surface marker 
expression to verify their phenotype (Figure 86). 

3D DNA FISH was performed as described previously [79] with 
BAC clones RP23-234A12 and RP23-435I4 (located at the distal 
end of the Vk region and at the Ck/ enhancer region, respectively; 
Figure 3A) obtained from BACPAC Resources (Oakland, CA). 
Probes were directiy labeled with Chromatide Alexa Fluor 488-5 
dUTP and Chromatide Alexa Fluor 568-5 dUTP (Invitrogen) 
using Nick Translation Mix (Roche Diagnostics GmbH). 

Cultured primary cells were fixed in 4% paraformaldehyde, and 
permeabilized in a PBS/0.1% Triton X-IOO/O.!""!) saponin 
solution and subjected to liquid nitrogen imm{;rsion following 
incubation in PBS with 20% glycerol. The nuclear membranes 
were permeabilized in PBS/0.5% Triton X-100/0.5% saponin 
prior to hybridization with the DNA probe cocktail. Coverslips 
were sealed and incubated for 48 h at 37°C, washed, and mounted 
on slides with 10 |J.l of Prolong gold anti-fade reagent (Invitrogen). 

Pictures were captured with a Leica SP5 confocal microscope 
(Leica Microsystems). Using a 63x lens (NA 1.4), we acquired 
images of ~70 serial optical sections spaced by 0.15 Jim. The 
datasets were deconvolved and analyzed with Huygens Profes- 
sional software (Scientific Volume Imaging, Hilversum, the 
Netherlands). The 3D coordinates of the center of mass of each 
probe were transferred to Microsoft Excel, and the distances 
separating each probe were calculated using the equation: ^(Xa— 
Xb)2+(Ya-Yb)2+(Za-Zb)2, where X, Y, and Z are tiie coordi- 
nates of object a or b. 

Statistical Analysis 

Statistical significance was analyzed using a nonparametric 
Mann-Whitney U test (IBM SPSS Statistics 20). The p values< 
0.05 were considered significant. 

Accession Numbers 

3C-seq and microarray expression datasets have been submitted 
to the Sequence Read Archive (SRA, accession number 
SRP032509) and Gene Expression Omnibus (GEO, accession 
number GSE53896), respectively. 

Supporting Information 

Figure SI Gene distance matrix analysis using gene 
expression profiling data from pre-B /pro-B cell frac- 
tions representing a pre-BCR signaling gradient. Micro- 
array expression profiling was performed on three or four 
independent FACS-purified B220"'CD19"' pre-B cell fractions 
from wild-type (WT), Btk, and Slp65 single- and double-deficient 
VH81x transgenic Ragl~'~ mice (see Figure 1 for gating strategy). 
The TIGR Multi Experiment Viewer software package (MeV 
version 4.8.1) was used to perform a one-way ANOVA analysis 
(/)<0.01) and identify genes differentially expressed within the five 
B-cell fractions (versus VH81X Tg Ragl~'~ pre-B cells). The 
software was subsequently used to create a gene distance matrix of 
highly significant genes, resulting in the depicted plot. Differences 
in gene expression profiles are depicted as a color code; darker 
colors indicate greater similarity, and brighter colors less similarity 
between groups. Consistent with the unsupervised clustering 
analysis shown in Figure IB, a clear gene expression gradient 
among the five B cell groups emerges in which Btk ' Slp65 ' 
pre-B and Ragl~'~ pro-B cells show highly comparably expression 
signatures. 
(TIF) 

Figure S2 Locus-wide 3C-Seq analysis of the IgVi 
region (Sis viewpoint) plotted as line graphs. Overview 
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of long-range interactions revealed by 3C-Seq experiments 
performed on the indicated cell fractions, representing a gradient 
of pre-BCR signaling. Shown are the relative interaction 
frequencies (average of two replicate experiments) for the Sis 
viewpoint per 100 kb region across the IgK locus, plotted as line 
graphs. The bottom graph shows an overlay of the WT and 
Btk ^ Slp65 ' pre-B cell interaction frequencies. The ~8.4 Mb 
region containing the Igfi locus (yellow shading) and flanking 
regions (cyan shading) is depicted. Pre-B cell fractions were FACS- 
purified from the indicated mice on a VH81x transgenic Ragl~'~ 
background (see Figure 1 for gating strategy). 
(TIF) 

Figure S3 Locus-wide 3C-Seq analysis of the Igvi region 
(IEk viewpoint) plotted as line graphs. Overview of long- 
range interactions revealed by 3C-Seq experiments performed on 
the indicated cell fractions, representing a gradient of pre-BCR 
signaling. Shown are the relative interaction frequencies (average 
of two replicate experiments) for the IEk viewpoint per 100 kb 
region across the Iga locus, plotted as line graphs. The bottom 
graph shows an overlay of the WT and !ilp65 ' pre-B cell 

interaction frequencies. The ~8.4 Mb region containing the Ign 
locus (yellow shading) and flanking regions (cyan shading) is 
depicted. Pre-B cell fractions were FACS-purified from the 
indicated mice on a VH81x transgenic Ragl~ ~ background (see 
Figure 1 for gating strategy). 
(TIF) 

Figure S4 Locus-wide 3C-Seq analysis of the /^K region 
(3 'Ek viewpoint) plotted as line graphs. Overview of long- 
range interactions revealed by 3C-Seq experiments performed on 
the indicated cell fractions, representing a gradient of pre-BCR 
signaling. Shown are the relative interaction frequencies (average 
of two replicate experiments) for the 3'Ek viewpoint per 100 kb 
region across the Ign locus, plotted as line graphs. The bottom 
graph shows an overlay of the WT and Btk~ ' ~ Slp65~ ' ~ pre-B cell 
interaction frequencies. The ~8.4 Mb region containing the Igti 
locus (yellow shading) and flanking regions (cyan shading) is 
depicted. Pre-B cell fractions were FACS-purified from the 
indicated mice on a VH81x transgenic Ragl background (see 
Figure 1 for gating strategy). 
(TIF) 

Figure S5 Selected zoom-in pictures of the 3C-seq data 
in the .^K locus and its upstream and downstream 
regions. (A) Map of the ~8.4 Mb genomic region containing the 
Iga locus (yellow shading) and flanking regions (cyan shading). 
Chromosomal coordinates and gene and viewpoint locations are 
also depicted. Locations of the zoom-in regions shown in (B) are 
represented as colored rectangles (red, upstream; green, down- 
stream; dark grey, Ign locus). (B) Average interaction frequencies 
per Bglll fragment (average of two replicate experiments) of the 
three regulatory elements with selected regions upstream, 
downstream, and within the Iga locus. Data from the five B cell 
precursor fractions representing a pre-BCR signaling gradient are 
shown. Complete locus-wide 3C-seq data can be found in Figure 3 
(plotted per individual Bglll fragment) or Figures S2, S3, and S4 
(plotted per 100 kb region), enh., enhancers. 
(TIF) 

Figure S6 Phenotypic characterization of cultured pre- 
pro-B, pro-B, and pre-B cells. Representative FACS analysis 
of CD43/CD19 (A) and CD2/CD19 (B) surface marker 

expression on IL-7 cultured E2a~'~ pre-pro-B, Ragl~'~ pro-B, 
and VH81x Ragl~'~ pre-B cells prior to 3D-FISH experiments. 
(TIF) 



Figure S7 Pre-BCR signaling is associated with an 
increase in the ratio of interactions inside the .^K locus 
over interactions outside the IgK. locus. The differential 

effects of pre-BCR signaling on long-range chromatin interac- 
tions of the IEk, 3'Ek, and Sis elements, as measured in 3C-seq 
datasets, were quantified (see Figure 4A). For aU three regulatory 
elements, the jv-axis shows the ratio of the average interaction 
frequencies per Bglll fragment inside the ~3.2 Mb IgK locus over 
the average interaction frequencies per Bglll fragment in the 
flanking regions (~2.0 Mb upstream together with 3.2 Mb 
downstream). Analyses of the five groups of B ceU precursors, 
representing a gradient of pre-BCR signaling, are shown, 
revealing that pre-BCR signaling is associated with a preference 
for interaction with fragments inside the Vk region over 
fragments outside the region. 
(TIF) 

Figure S8 Interaction frequencies of K regulatory ele- 
ments with nonfunctional V,^ genes as measured in B cell 
fractions representing a pre-BCR signaling gradient. 

Quantitative analysis of 3C-Seq datasets obtained for the five B 
cell precursor fractions representing a pre-BCR signaling gradient, 
using the three indicated K regulatory elements as viewpoints. 
Average interaction frequencies within the Vk region were 
determined for fragments that do not contain any gene ("no 
Vf^"), those that contain a nonfunctional Vk gene {"pseudo V^"], and 
those that contain a functional Vk gene {"Junctional Vk"). See 
Materials and Methods section for more details on analysis 
methods. Statistical signiffcance was determined using a Mann- 
Whitney U test (n.s., not significant, ^S0.05). 
(TIF) 

Figure S9 Correlations between the interaction 
profiles of the three K regulatory elements. Correlation 
plots of average interaction frequencies of the three regulatory 
elements with the 101 functional Vk genes (A) or only the Vk3 
gene family (B) are shown for WT pre-B cells (top, gray labels) 
versus Btk ' Slp65 ' pre-B cells (bottom, yellow labels). Note 
that under conditions of low pre-BCR signaling {Btk~ ' ~ Slp65~ ' ~ 
pre-B cells) correlation strength is significantiy reduced, with the 
exception of the correlation between de Vk3 interaction profiles of 
the Sis and IEk elements. 
(TIF) 

Table SI Genes down-regulated in the absence of Btk 
and Slp65. 

(DOC) 

Table S2 Genes up-regulated in the absence of Btk and 
Slp65. 

(DOC) 

Table S3 Binding of E2a and Ikaros to the K enhancers 
and V,^ genes in wild-type and Btk or Slp65-deficient pre- 
B cells. 

(DOC) 
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